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DELAY AND CHANNEL ESTIMATION FOR 
MULTI-CARRIER CDMA SYSTEM 

BACKGROUND OF THE INVENTION 

This invention relates to receivers for decoding signals received via multiple 
propagation paths having different propagation delays and, more particularly, to delay and channel 
estimation for multi-carrier systems. 

Radio receivers are often used for decoding fading signals with the aid of estimates 
of the instantaneous propagation channel phase and amplitude characteristics. An exemplary 
application for such a radio receiver is a cellular phone for decoding signals transmitted by base 
stations using code division multiple access (CDMA) protocols. 

A radio receiver, such as for a CDMA system, receives digitally coded and modulated 
signals from a transmitter. These signals include known, preselect signal patterns at known time 
intervals. Using known signal patterns, also referred to as pilot sequences, or pilot channel as 
commonly used in the CDMA literature, the receiver forms successive estimates of the delay and 
amplitude or complex value of propagation path characteristics between the transmitter and the 
receiver. These include estimates for multiple paths in the case of multi-path propagation. 

In a typical CDMA system a signal is transmitted on a single carrier. However, this 
can limit data rate of transmission. More recently, a multi-carrier (MC) system is proposed in which 
the signal is transmitted using three sub-carriers. This effectively triples the data rate. However, the 
MC system requires that filtering , correlation and de-modulation be performed for each of the three 
sub-carriers. 

The common approach for MC CDMA system delay and channel estimation is to use 
three baseband filters. The signal for the first sub-carrier is extracted using the first filter. A first 
correlator correlates the filtered signal using the known sequence. From this, the delay and complex 
gain can be determined. This process is repeated for the second and third sub-carriers. 
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In general, the signal transmitted by the k'th carrier is given by 



(i) 



where s/t) is the known pilot signal of the k'th carrier, o/t) is the sum of all other Walsh-orthogonal 
channels in the same carrier and e^ h is a constant phase term. Both s^i) and o^t) are baseband 
signals resulting from modulating data sequences on to streams of baseband pulse shape. The overall 
signal sent by the transmitter is 

«(o=Z«*('y 2 * 04+/,) ' (2) 



where f„ is the center carrier frequency and f t is the frequency separation of the fc'th carrier with 
respect to the center carrier frequency. The frequency spacing^' s are assumed wide enough so that 
there is no spectrum aliasing between sub-carriers. 

After propagating through a multi-path channel, the RP received signal can be 

expressed as 

? r (0=Zc / g(r-r / ) + «(f) 

(3) 



-1c, 



+ »('), 



where c 7 is the complex channel gain of the multi-path, r, is its delay and n(t ) is the Additive White 
Gaussian Noise (AWGN) with power spectral density A^. 
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SUMMARY OF THE INVENTION 

In accordance with the invention a multi-carrier CDMA receiver identifies multi-paths 
and relative delays by considering the multi-carrier signal as an entire wide band signal and performs 
estimation using the known signals in all sub-carriers. 

Particularly, The receiver receives a signal transmitted on plural sub-carriers each 
having a known pilot sequence>A,plurality of down-converters down-convert the received signal to 
baseband signals. A delay and channel>stirnator correlates the baseband signals to produce an 
estimate of channel gain and multi-path delay. A plurality of demodulators, one for each of the plural 
sub-carriers, is coupled to the delay and channel estimator, eacKjiemodulating one of the baseband 
signals using the estimate of channel gain and multi-path delay. 



BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram of a mobile terminal including a receiver according to the 

invention; 

Fig. 2 is a block diagram of a prior art receiver; 

Fig. 3 is a curve illustrating auto correlation functions for the prior art receiver and 
a receiver according to the invention; 

Fig. 4 is a block diagram of a receiver according to one illustrated embodiment of the 

invention; 

Fig. 5 is a block diagram of a receiver according to another illustrated embodiment 
of the invention; 

Fig. 6 is a set of graphs illustrating a method of synthesizing a radio channel profile 
in accordance with the invention; and 

Fig. 7 is a block diagram of a mobile communication system including a receiver 
according to the invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

Fig. 1 is a block diagram of a mobile terminal shown generally at 10. The mobile 
terminal 10 includes an antenna 12, a receiver 16, a transmitter 18, a speaker 20, a processor 22, a 
memory 24 a user interface 26 and a microphone 32. The antenna 12 is configured to send and 
receive radio signals between the mobile terminal 10 and a wireless network (not shown). The 
antenna 12 is connected to a duplex filter 14 which enables the receiver 16 and the transmitter 18 to 
receive and broadcast (respectively) on the same antenna 12. The receiver 16 demodulates, 
demultiplexes and decodes the radio signals into one or more channels. Such channels include a 
control channel and a traffic channel for speech or data. The speech or data are delivered to the 
speaker 20 (or other output device, such as a modem or fax connector). 

The receiver 16 delivers messages from the control channel to the processor 22. The 
processor 22 controls and coordinates the functioning of the mobile terminal 10 responsive to 
messages on the control channel using programs and data stored in the memory 24, so that the mobile 
terminal 10 can operate within the wireless network. The processor 22 also controls the operation 
of the mobile terminal 10 responsive to input from the user interface 26. The user interface 26 
includes a keypad 28 as a user-input device and a display 30 to give the user information. Other 
devices are frequently included in the user interface 26, such as lights and special purpose buttons. 
The processor 22 controls the operations of the transmitter 18 and the receiver 16 over control lines 
34 and 36, respectively, responsive to control messages and user input. 

The microphone 32 (or other data input device) receives speech signal input and 
converts the input into analog electrical signals. The analog electrical signals are delivered to the 
transmitter 18. The transmitter 18 converts the analog electrical signals into digital data, encodes the 
data with error detection and correction information and multiplexes this data with control messages 
from the processor 22. The transmitter 18 modulates this combined data stream and broadcasts the 
resultant radio signals to the wireless network through the duplex filter 14 and the antenna 12. 

The mobile terminal 10 may be used in a multi-carrier (MC), code division multiple 
access (CDMA) system in which the signal is transmitted using plural, such as, for example, three 
sub-carriers. Particularly, the receiver 16, in conjunction with the processor 22 utilizes the multi- 
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carrier signal as an entire wideband signal and performs channel estimation using the known pilot 
signals in all sub-carriers. 

The conventional approach for multi-carrier CDMA system delay and channel 
estimation is solved using the single-carrier method, as discussed above. Fig. 2 is a block diagram 
of such a conventional system. Equivalently, for the i*th sub-carrier the RF signal received on the 

antenna 12 is first down-converted to baseband by multiplying it withe~ j/2,r( ^ 0+ ^ ) ' using multipliers 
40-/. For simplicity hereinafter, elements common to each sub-carrier are labeled in the drawings with 
a suffix -/, where /=1,2 or 3, corresponding to the respective carrier, but referred to herein with the 
suffix omitted. The output of each multiplier 40 is filtered using a baseband filter 42 to remove the 
other sub-carriers. Then delay and channel estimation are performed independently on each sub- 
carrier using the single-carrier method with correlators 44. The correlators 44 estimate the channel 
gain and multi-path delay. This information, along with the filtered signal, is applied to demodulators 
46. The demodulators 46 use the delay and channel estimates to demodulate the received signal for 
each sub-carrier in a conventional manner. Demodulated data is then sent to a decoder (not shown). 

After down conversion, the resulting baseband signal for the /*th sub-carrier is 

= X cfi^^e^^u.it - r,) (4) 



k*i 



where nft) is the down-converted Additive White Gaussian Noise (AWGN). Estimation of c l and r, 

can be achieved by correlating the received signal rfi) with a replica of the /'th known pilot signal 
with a hypothesized delay t. Since the pilot signal is a baseband signal, the correlation filters out the 

signals of other carriers, i.e., the second term in Eq. (4). Replacing c^ 2 ^ 1 with c { for 
simplification, the output of the correlation as a function of x can then be expressed as 
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liT)=\r i (t)s l {t-T)e»'dt 

= \ {Z c,e J1 ' f '\ S ,{t- r,)+ o,(t- T,)]e»' + «,(') 
= £ cfi-»""\R£T- r,)+ /^(r- r,)] + »,(r), 



s]{t- r)e' J *'dt 



(5) 



where 

i? J( (r)= Jj,(Oj;(f-r>* 

is the auto-correlation function of sfi), 

is the cross-correlation between o/y and j/y and 
is the filtered noise. 

By varying the value of t within the observation period, the delay profile of the 
channel as a function of t can be obtained. The delays with correlation output exceeding a certain 
threshold are identified as possible multi-path echos. The channel gain of the echo is simply its 
corresponding correlation output. In particular, the correlation output of the /w'th echo is 

^(0= V y2,r/ ' r "K(0)+^, Oi (0)] 

+ I ^""'k (r. - r f ) +R„ (r m - r,)] + n,(rj (6) 

l*m L J 
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where E h = R Sj (0) = J |s ( (t)\ dt is the energy of the pilot signal and 

= I *i'" yJ "' r 'k('« ~ ',) + " r,)l (7) 

Urn L J 

is the interference from other multi-paths. The cross-correlation term R s (0) is removed due to 

the orthogonality between users in the same sub-carrier. In a typical CDMA system, the multi-path 
interference can be approximated by AWGN. Therefore, the performance of the channel estimation 
depends on the Signal to Noise Ratio (SNR) between the pilot signal and the sum of the multi-path 
interference and AWGN. The performance of the delay estimation, on the other hand, depends on 
the main lobe of the auto-correlation R s ( r ) of the pilot signal. The auto-correlation function with 

narrower main lobe has higher time resolution and therefore better performance in delay estimation. 
If the interference and noise terms in Eq. (6) are small, 

A ,('.}' » c m e-'*""- , (8) 

which gives the channel estimate of the /w'th path in the /'th carrier. The effect of fading can then be 
reversed by multiplying the received signal with the complex conjugate of Eq. (8) and the coherent 
demodulation of the information bits in oft) can be achieved. 

The first step of the system according to the present invention maybe to down-convert 

the received signal in Eq. (3) to baseband (with respect to f 0 ) by multiplying it with e' J2n ^ 0 ' . The 
resulting baseband signal is 

r(/)=r(/)g- J2 * /i ' 

(9) 

k 



From Eq. (9) it is clear that to fully utilize the entire wideband signal, the known signal to correlate 
the receive signal r(t) with should be 
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(10) 



Thus, the correlation output as a function of z is 



dt 



k 

= £ e n * fkX \ [r(t)e- J2ff ^sl(t- r)e jh ]dt 

k 



00 



Substituting Eq. (5) in Eq. (1 1) gives 
* (0= I ^""{l C ,<r^"[^(r- r,) + /^(r- r,)]+ n k {r)> 

= I c,S {^"'"""'[^(r- r,) + ^, t (r- r,)]} (12) 

Comparing Eq. (12) with its single carrier counterpart in Eq. (5), it can be noted that 
the auto-correlation function dictating the estimation performance is now 

^,(o=Z« y2 * /i, ^(o. 03) 



Similarly, the correlation output of the m'th echo is 
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(14) 



It is clear from Eq. (14) that if there are K sub-carriers in the system and all sub- 



carriers have the same power, the disclosed approach has a 10 \og(K) dB SNR gain over the 
conventional approach discussed above relative to Fig. 2. Furthermore, the accuracy of the delay 
estimation depends on the new composite auto-correlation function given in Eq. (13), which has a 
main lobe that is \/K of that of a single carrier since the composite pilot signal has K times the 
bandwidth of the single carrier (A W. Rihaczek, "Principles of High-Resolution Radar," Artech 
House Radar Library, 1996). 



where q K * } is a random sequence of QPSK (quadrature phase-shift keying) symbols, p(t) is a band- 
limited pulse- shaping signal generally referred to as a chip, and T c is the chip interval. The pilot 
channels in theMC mode of IS-2000 (TIA/EIA/IS-2000.2 [Ballot Version], "Physical Layer Standard 
for cdma2000 Spread Spectrum Systems," 1989) can be effectively modeled by Eq. (15) with \/T c 
- 1.2288 MHz and an approximate Nyquist pulse shaping function. There are a total of three sub- 
carriers separated by 1.25 MHz in the system. Fig. 3 shows the real part of the auto-correlation for 
one of the single carrier pilot signals (IX) and that of the multi-carrier composite pilot signal (3X). 
A section of 128 chip duration is extracted randomly as the pilot signal and the SNR defined as the 
ratio of the energy per chip per carrier over N 0 is set to -7 dB. The resolution improvement in delay 
estimate is apparent from the narrower peak width of the 3X function which is centered at zero delay, 
and SNR improvement for channel estimate is evident from the higher peak of the 3X function. 



As an example, consider the pilot signal 




(15) 



ft 



It is clear from Eq. (11), which is repeated here without the intermediate derivation 



steps, 




(16) 



k 
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that the correlation output according to the present invention can be obtained from the individual 
correlation outputs of the conventional approach of Fig. 2. In other words, a channel profile with 
higher resolution and stronger SNR can be synthesized from several of those with lower resolution 
and weaker SNR using Eq. (16). 

The channel estimate derived from Eq. (16) at the delay x m is c m . For the coherent 

demodulation of the j'th carrier, this coefficient needs to be rotated by a factor of e~ i2n f,Tm . In cases 
where the coherent bandwidth of the fading channel is of the order of the bandwidth of a single 
carrier, the high resolution correlation output X(r) will likely see more than one path in a chip period 
while the individual correlation output A,(r) will see only a blurred image of one single path resulting 
from the superposition of several closely separated paths. The channel coefficients as observed by 
each individual sub-carrier will therefore differ not only in phase, but also in amplitude. However, 
they can still be derived from Eq. (16) with appropriate interpolation if needed. 

In cases where the coherent bandwidth of the channel is so small that the interpolation 
of the higher-resolution channel profile into the lower-resolution channel profiles does not improve 
the performance of the individual sub-carriers or does not justify the computational complexity, a 
switch function can be implemented in the processor to decide whether to use the conventional 
approach or the disclosed approach based on the measurement of the channel's coherent bandwidth. 

The disclosed approach according to the present invention identifies all of the pilot 
(known) signals in the sub-carriers as a single wideband pilot signal and correlates it against the entire 
multi-carrier received signal. The resulting correlation output outlines the profile of the multi-path 
channel with higher temporal resolution and stronger SNR. 

There may be several embodiments of the system according to the present invention. 
A receiver architecture implementing the above general description is shown in block diagram form 
in Fig. 4. Where elements are similar to those shown in Fig. 2, like reference numerals are used. 
Initially, the received signal is converted to baseband with multipliers 40 and filtered using filters 42, 
as with the conventional approach of Fig. 2 The difference from the conventional approach is 
summarized in the following: 
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The RF signal is down-converted to baseband with respect to the center carrier 
frequency^ by multiplying the received signal withe' 2 *^ 0 ' using a composite multiplier 50: 

/•(/) = r(t)e' j2fffot . (17) 

The composite baseband signal is filtered with a filter 51 that passes all of the sub-carriers. Channel 
estimation and delay are determined by a correlator 52. The correlator 52 correlates the down- 
converted wideband signal with 5* (f) = (t)e~ j2jI ^ e~ J ** as follows: 

X ( T )= jr(t)s*(t- r)dt (18) 

From k(r), a multi-path can be identified and its relative delay f m estimated using a conventional 
technique such as threshold comparison. The channel coefficient estimate of this path is simply 

C m = (19) 

Finally, C m is rotated bye~ jlnftTm to give the channel estimate for the /'th carrier at 

delay x m . The channel delay and estimates, along with the sub-carrier baseband signals from the filters 
42, are then supplied to the demodulators 46, as with the conventional approach. 

Eq. (16) suggests an alternative embodiment of a receiver architecture illustrated in 
Fig. 5. It differs from the first embodiment of Fig. 4, only in the calculation of A(t). 

As with the conventional approach, the alternative embodiment down-converts the 
RF signal of each individual sub-carrier to baseband with respect to its own center frequency by 

multiplying the received signal with e~ j2iT ^ 0+ ^' above using multipliers 40, as follows: 

ri {t)^r{t)e' j2n{f ^ f,) \ (20) 



The signals are then filtered with filters 42 to provide the sub-carrier baseband signals. For the fth 
sub-carrier, a correlator 44 correlates the down-converted signal with s*{t)e~ j ^ : 
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*,(0= Jr^fKCr-r)^**. (21) 

In accordance with the present invention, the high resolution channel profile A(t) is 
synthesized by combining the individual correlation outputs. This is done by rotating the individual 
correlator outputs using multipliers 54. The signals are then summed and derotated at a block 56. 
This combination is expressed as 

,1 (r)=5> 2 ^,(0 (22) 

From X(t), a multi-path can be identified and its relative delay i m estimated using a conventional 
technique such as threshold comparison. The channel coefficient estimate of this path is simply 

c m = ^ V mJ . (23) 

Finally, c m is derotated by e~ J2 * ^ ,Tm in the block 56 to give the channel estimate for the fth carrier 
at delay x m . 

As is conventional, the filter 42 is implemented in hardware in the receiver 1 6, see Fig. 
1. The equations described herein, and represented by blocks in the block diagrams of Figs. 4 and 
5 are implemented in software by the processor 22. However, the invention is not limited to such a 
division of functionality. As such, where the term receiver is used herein, the term encompasses 
functions implemented in hardware in the receiver block 16 and in software in the processor block 
22. 

The embodiment of Fig. 4 uses all of the known signals in the sub-carriers as a single 
wideband reference signal and correlates this reference signal against the entire multi-carrier received 
signal. The resulting correlation output yields the underlying channel profile with stronger SNR and 
higher temporal resolution. Furthermore, this higher resolution channel profile can be synthesized 
by combining the individual lower resolution profiles of the sub-carriers using equation (11). 
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While equation (1 1) is expressed in continuous time, the correlations are performed 
on discrete samples in practice. Therefore, if the individual correlation output X k (x) is sampled every 
t a seconds, then a straightforward method for synthesizing and sampling A(t) is 

k 

which has the same sampling rate as the individual sampled correlation output X(nt,). 

In accordance with another aspect of the invention a receiver solves the problem of 
sampling and synthesizing the composite correlation output, or k(x) as expressed in mathematical 
term described above, required for improving channel estimation in a multi-carrier system. As will 
be described, the receiver can synthesize X(t) at any desired sampling rate if the sampling rate of the 
individual correlation outputs meets the Nyquist criterion. 

Although very simple, equation (24) is not an optimal way of combining several time 
sequences from different frequency bands. In the frequency domain, the Fourier transform of A(t), 
denoted by A(f), is simply 

A [f) k \x(r)e^'dT =lk k (f~ /.) (25) 

where \(J) is the Fourier transform of A k (t) centered at zero frequency. When sampled at or above 
the Nyquist rate (twice the bandwidth of A^if)), A k (nt,) contains sufficient information to reconstruct 
the continuous time signal A. k (x), and thus also its frequency-offset version e"^ r A k (t). Therefore, the 
composite correlation output A(x) can be reconstructed fully and sampled at any desired rate if X k (x) 
is sampled at or above Nyquist rate for all sub-carriers. 

Furthermore, in order to reduce computational complexity, X k (t) is usually sampled 
at the Nyquist rate. Thus, equation (24) will not meet the Nyquist criterion since X(x) has K times 
the bandwidth of the individual spectrum for a system with K sub-carriers. To satisfy the Nyquist 

K 

criterion, A(t) must be sampled at or above the rate of — . The following explains how to 
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synthesize A(x) at the desired sampling rate using Discrete Fourier Transform (DFT) from X k (x)'s, 
which are sampled at or above the Nyquist rate. The sampling may be done, for example, at two 
samples per chip or greater. 

This aspect of the invention is explained by an example as shown in Fig. 6. This 
example relates to the receiver architecture of Fig. 5. Let 



denote the discrete samples over a period of Nv s of the individual correlation output sampled at the 
rate of 1/t b Hz, which is assumed to meet the Nyquist criterion. Graphs 60, 62 and 64 represent the 
individual correlation outputs for respective first, second and third sub-carriers in a three carrier 



the same observation period by the following procedure: 

1 . For each sub-carrier k = 0, 1 , . . K - 1 , compute the DFT of length N, denoted by 



A k [n]=A k (nT,)> for n = 0, 1, N- 1 



(26) 



L 



system. The composite correlation output A(x) sampled at the rate of 



Nr. 



Hz can be obtained for 



F„{.},forA k [ W ]: 





(27) 



over the frequency range of m - -A/72, -Nil + 1, 0, 1, Nil • 1. The results are 
represented by graphs 66, 68 and 70. 



2. Compute m h the carrier frequency offset f k in discrete domain, for each sub- 
carrier: 

m k =f k Nr s . (28) 
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3. Form the discrete Fourier transform of X(z) by summing all the individual DFT 
spectra, shown in graphs 66, 68 and 70, shifted by the appropriate frequency offset: 

AM=Z A (29) 



for m - -L/2 t -L/2 + 1, 0, 1, L/2 - 1. Zeroes are added if necessary to give a 
total of L samples. The resultant transform is represented by the graph 72. 



4. Calculate the inverse DFT of length L for A[m]: 

^%FT{ A M} (30) 

L/2— 1 jljtmn 
m=-L/2 

The resulting time domain sequence X[n], represented by the graph 74, is the 
composite correlation output sampled at the desired rate. 



As is apparent, the resulting time domain sequence has a higher resolution profile 
illustrated by the graph 74 having two peaks at substantially higher correlation amounts compared 
to the individual correlation output curves 60, 62 and 64. 

The receiver embodiments discussed above relate to the mobile terminal 10. As is 
apparent similar receiver embodiments could be used in a receiver 76 of a base station 78, as shown 
in Fig. 7. The base station 78 communicates with the mobile terminal 10 via a mobile communication 
network, represented at 80. The receiver 76 would be used in applications where the mobile terminal 
10 transmits a multi-carrier signal. 



